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I present the results of first principles calculations of the phonon dispersions and electron-phonon 
coupling for BaTi2Sb20. The phonon dispersions show a weak lattice instability near the zone 
corners that leads to a CDW phase. The calculations of the electron-phonon coupling reveal strong 
coupling, especially to the in-plane Ti modes. The total coupling is large enough to readily explain 
the superconductivity in this compound. As the Fermi surfaces are disconnected with different 
orbital character weights, this compound is likely to host a multiband superconductivity. The role 
of spin fluctuations in suppressing the superconducting Tc is also discussed. 

PACS numbers: 74.25. Kc,74.20.Pq,74.70.-b 



C/3 . 



I. INTRODUCTION 

The discovery of high temperature superconductivity 
in iron based compounds by Kamihara and coworkers^ 
has increased interest and intensified efforts to find super- 
conductors in new families of compounds. These efforts 
have focused on searching for materials with structural, 
electronic and magnetic properties similar to that of the 
iron or cuprate superconductors. In particular, the ideal 
candidates are thought to be materials with a layered 
structure that are in proximity to a spin density wave 
instability due to Fermi surface nesting like in the iron 
based superconductors or a Mott insulating phase due to 
strong on-site Coulomb repulsion as in the cuprates. 

Recently, Doan et al. have reported superconductivity 
in Bai_2;Na2;Ti2Sb20 with a maximum Tc of 5.5 K.— Su- 
perconductivity in this compound has also been reported 
by Yajima et al. with a Tc of 1.2 K4 In addition to the 
superconducting transition, there is also an anomaly at 
= 54 K that manifests by showing strong features in 
the measurements of susceptibility, resistivity, and spe- 
cific heat. The microscopic mechanism and the order 
parameter for this transition has not been determined, 
but it has been ascribed to cither a charge density wave 
(CDW) or spin density wave (SDW) instability. Fur- 
thermore, the superconducting phase coexists with this 
CDW/SDW phase. 

This compound occurs in a tetragonal structure with 
spaccgroup PA/mmm and consists of square planes of 
transition element Ti. The O atoms are placed alter- 
nately at the center of the Ti squares such that each 
Ti atom has two O nearest neighbors. The Sb atoms 
are placed above and below the center of the Ti squares 
that do not contain O atoms. The Ti2Sb20 layers are 
stacked alternately with layers of Ba atoms along the c 
axis. This layered structure is similar to the iron based 
or cuprate superconductors as all of them have square 
planes of transition metal atom as a structural motif. 

In addition, Singh has further highlighted similari- 
ties in the electronic and magnetic properties between 
BaTi2Sb20 and the iron superconductors based on first 



principles calculations^^ The Fermi surface shows sub- 
stantial nesting and leads to an antiferromagnetic insta- 
bility as in the case of iron superconductors. This nested 
Fermi surface is consistent with the results of Pickett on 
a similar compound Na2Sb2Ti20, although he was un- 
able to find a magnetically ordered statci^ The electronic 
structure of Na2Ti2Sb20 was also studied by de Biani 
et al. using tight-binding calculations, and they also find 
substantial nesting in the Fermi surface that may lead to 
a CDW instability,^ 

In any case, the presence of a nested Fermi surface and 
antiferromagnetic instability in BaTi2Sb2 leads Singh to 
predict a spin fiuctuation mediated superconductivity 
similar to that of the iron based superconductors ,i The 
pairing state is predicted to have a sign-changing s-wave 
symmetry that is different from that of the iron based 
superconductors. Furthermore, spin fluctuations are re- 
pulsive in the singlet channel and hence are strongly pair- 
breaking for the attractive electron-phonon interaction. 
Thus, it is reasonable to expect that the superconductiv- 
ity in BaTi2Sb20 is mediated by spin fluctuations, espe- 
cially if the phase transition at Ta = 54 K that coexists 
with the superconducting phase is due SDW instability. 

However, there are some differences between 
BaTi2Sb20 and other superconductors near mag- 
netism. The maximum Tc — 5.5 K of this compound 
is significantly smaller than that of the iron based or 
cuprate superconductors. This could be due to a small 
value of Stoner parameter for Ti^ or pair-breaking effects 
of electron-phonon interaction. It is worthwhile to note 
that Sr2Ru04, which is purported to be an unconven- 
tional superconductor;^ also has a small Tc of 1 KJ^ A 
more glaring difference is the very sensitive dependence 
of superconductivity on other constituent elements 
of the compound. A transition similar to the one in 
BaTi2Sb20 is also seen in other members of the family 
that includes Na2Ti2Pn20 {Pn = As, Sb), BaTi2As20, 
(SrF)2Ti2P7i20 {Pn = As, Sb), (SmO)2Ti2Sb20^-iS 
but superconductivity has so far only been observed 
in Na doped BaTi2Sb20. This is in contrast to iron 
superconductors which exhibit superconductivity for a 
wide variety of cations and pnictogens. In the cuprates. 
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superconductivity is also found for a wide variety of 
interlayer fillings that separate the Cu02 layers. Even in 
the ruthenates, doping at the cation site or changing the 
interlayer spacing reveals various magnetic interactions. 
This suggests that magnetism and mechanism for super- 
conductivity in doped BaTi2Sb2 is quite different from 
other unconventional superconductors, and perhaps 
lattice instabilities and clcctron-phonon coupling play 
an important role. 

In this paper, I present the results of first principles 
calculations that show presence of a CDW instability 
and a total electron-phonon coupling strong enough to 
yield a conventional superconductivity in BaTi2Sb20. 
The phonon dispersions reveal a weak lattice instabil- 
ity near the Brillouin zone corners. This instability is 
associated with the elongation or compression of the Ti 
squares without an enclosed O such that the Ti squares 
with an O rotate cither clockwise or counterclockwise by 
a small amount. I also find presence of strong clcctron- 
phonon couplings, especially near the zone corners. The 
total electron-phonon coupling Aop = 1.28 is large enough 
to readily explain the superconductivity in this material. 



II. METHODS 

The phonon dispersions and clcctron-phonon results 
presented here were obtained using density functional 
perturbation theoryi^ within the generalized gradient ap- 
proximation of Perdew, Burke, and Ernzerho^ as imple- 
mented in the Quantum-ESSPRESSO package.— I used 
the experimental lattice parameters (a = 4.11039 A, c 
= 8.0864 A)^ but relaxed the internal Sb height param- 
eter zsb. I obtained a value for zsh = 0.2493, which 
agrees well with both the experimcntal^j^ and calculated^ 
values. This is in contrast the underestimation from 
density functional calculations found in the iron based 
superconductors.— This indicates that magnetism and/or 
lattice dynamics in BaTi2Sb2 might be different from 
that of the iron based superconductors. 

I used pseudopotentials that were generated with the 
electronic configurations 5s^5p^5dP6s^6p^ , 3s^3p^As^3d^ , 
5s'^5p^5d", and 2s'^2p'^3d-^ , respectively, for Ba, Ti, Sb, 
and O. In particular, it was necessary to include the semi- 
core states to get the electronic structure in close agree- 
ment with the results of the full-potential calculations. I 
used the cut-offs of 60 Ry and 600 Ry for basis set and 
charge density expansions, respectively. An 8 x 8 x 4 k- 
grid was used for the Brillouin zone integration during 
the self-consistency. The dynamical matrices were also 
calculated on an 8 x 8 x 4 q-grid, while the double-delta in- 
tegration in the calculation of the electron-phonon spec- 
tral function was done on a 24 x 24 x 12 grid. I used 
the generalized full-potential method as implemented in 
the WIEN2k package^^ to perform the relaxation of the 
internal coordinates of a {\/2 x x 1) superccU. The 
muffin-tin radii of 2.4, 1.7, 2.4, and 1.7 Bohr were used 
for Ba, Ti, Sb, and O, respectively, and an 8 x 8 x 8 k-grid 
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FIG. 1: Calculated GGA band structure of non-spin-polarized 
BaTi2Sb2 in the underdistorted structure with experimental 
lattice parameters but relaxed internal coordinates. 

was used for the Brillouin zone integration. 



III. ELECTRONIC STRUCTURE 

The electronic structure and magnetism of BaTi2Sb20 
has been well described by Singhi^ I summarize his re- 
sults here for the sake of completeness. The band struc- 
ture of this compound calculated using the pseudopoten- 
tials and plane-wave basis set is shown in Fig. [T] This 
agrees reasonably well with Singh's full-potential band 
structure that also includes the spin-orbit interaction. I 
also did full-potential band structure calculations with- 
out the spin-orbit coupling, and I find excellent agree- 
ment with the results from psuedopotential calculation. 
This shows that the pseudopotentials that I used are ro- 
bust. 

The three bands that lie between —7.3 and —4.7 eV 
have a dominant O 2p character, and another six bands 
with mostly Sb 5p character lie between —4.7 and —0.5 
eV, relative to the Fermi energy. The bands near the 
Fermi level have mostly Ti Sd character with some hy- 
bridization with Sb 5p states, while the Ba 6s states are 
above the Fermi level. This suggest that Ti ions are nom- 
inally trivalcnt and are in d^ state, although the actual 
electron count will be different due to some covalcncy 
with O 2p and Sb 5p states. 

The Fermi surface (not shown; see Ref.i) has a very 
two dimensional, square electron sheet around zone cor- 
ner M with mixed Ti d^^, d^^-y^, and dxy characters. 
There is a three dimensional hole sheet around X with 
mostly Ti d^^ with some admixture of c?^2_j,2 and dxy 
characters. Around the zone center there is a complex 
three dimensional electron sheet with mixed dz2 and dxy 
character. In the calculations without spin-orbit cou- 
pling, there is also a small electron sheet around the zone 
center that has a mixed character. Singh calculated the 
real part of the susceptibility within the constant matrix 
element approximation and finds peaks in the susceptibil- 
ity due to substantial nesting between Fermi sheets, espe- 
cially near the q vector (5,0) (27r/a). In addition, there is 
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FIG. 2: Calculated GGA phonon dispersions of non-spin- 
polarized BaTi2Sb2 in the underdistorted structure with 
experimental lattice parameters but relaxed internal coordi- 
nates. The imaginary frequencies are shown as negative. 



also a smaller peak near the q vector (^,5) (27r/a). This 
could lead to magnetic instabilities at these q vectors, al- 
though the height of the peaks could be diminished when 
the matrix element is explicitly taken into account in the 
calculation of the susceptibility. In any case, Singh finds 
that there is a magnetic instability at X that leads to a 
double stripe antiferromagnetic order with a doubling of 
the unit cell along either a or 6 axis with small moments 
of 0.2 /is per Ti. Based on the proximity to the SDW 
instability, Singh predicts a spin fluctuation mediated su- 
perconductivity in doped BaTi2Sb20 with a pairing state 
that has a sign-changing s-wave symmetry that is differ- 
ent from the one in the iron based superconductors. 



IV. PHONONS AND ELECTRON-PHONON 
COUPLING 

Although there seems to be a magnetic instability in 
BaTi2Sb20, this does not preclude a competing lattice 
instability or an electron-phonon superconductivity. In 
this section, I show results that indicate that the anomaly 
at Ta ~ 54 K is due to a CDW instability and the 
electron-phonon coupling is large enough to readily ex- 
plain the observed superconductivity. 

The calculated phonon dispersions of BaTi2Sb20 in 
the undistorted tetragonal structure with the experimen- 
tal lattice parameters and relaxed internal coordinates 
are shown in Fig. [2l and the corresponding phonon den- 
sity of states along with atomwise projections are shown 
in Fig. [3l A conspicuous feature of the phonon disper- 
sions is the presence of an unstable mode around M 
(^,5,0) and A indicating a CDW instability. 

The instability is weak with the maximum imaginary fre- 
quency of 30i cm~^, suggesting a very shallow double 
well potential. The unstable mode has only in-plane Ti 
character, and the distortion corresponds to elongation 
or compression of the Ti squares without an enclosed O 
such that the Ti squares with O rotate either clockwise 
or counterclockwise by a small amount (see Fig. 2]). I 
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FIG. 3: (Color online) Top: Calculated GGA phonon density 
of states G{uj) along with atomwise projections on Ba, Sb, and 
O of non-spin-polarized BaTi2Sb2 in the undistorted struc- 
ture with experimental lattice parameters but relaxed internal 
coordinates. Bottom: Phonon density of states weighted by 
projections on Ti atom as well as its in-plane and out-of-plane 
characters. 



also performed relaxation of a (-\/2 x \/2 x 1) super-cell 
with the full potential to check if the instability is a spu- 
rious effect of using pseudopotentials, and I again found 
a presence of this instability. The Ti ions get displaced 
by only 0.096 A from their original high-symmetry po- 
sition. Although the Ti displacements are very small, I 
find a substantial reduction of the electronic density of 
states near the Fermi level with a reduction from a value 
of 3.88 eV~^ per formula unit in the undistorted case 
to a value of 2.83 eV~^ per formula unit in the CDW 
phase, as obtained from full-potential calculations with- 
out spin-orbit coupling. In addition, I also performed 
frozen phonon calculations on the (a/2 x \/2 x 1) super- 
cell by displacing Ti atoms according to the eigenvector 
of the unstable phonon mode. As shown in Fig. [Sj the 
potential is anharmonic with a minimum at 0.082 A. and 
the depth of the well is very shallow 14.5 cm~^. Hence, 
it is not surprising that this CDW instability has so far 
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FIG. 4: (Color online) The distorted Ti plane in the CDW 
phase. The black grid shows the (\/2 x %/2 x 1) supercell 
relative to the undistorted structure. The big (cyan) balls 
represent Ti and small (red) balls are O. The Ti displacements 
are exaggerated to make the distortions conspicuous. 




FIG. 6: Calculated Eliashberg spectral function of 
BaTi2Sb2 in the underdistorted structure. The spectral 
weights for imaginary phonon frequencies are set to zero. 
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FIG. 5: Calculated anharmonic double-welled potential for 
the unstable mode obtained using frozen phonon method on 
a {^/2 X y/2 X 1) supercell. The energy is given per formula 
unit. 



been eluding detection. 

Returning back to the phonon dispersions, the 18 
modes of BaTi2Sb20 in the undistorted structure ex- 
tend up to 510 cm~^. There is a manifold of 16 bands 
with mixed Ba, Ti, Sb, and out-of-plane O character ex- 
tending up to 370 cm~^, which is separated by a gap 
from two in-plane O modes between 445 and 510 cm~^. 
Within the lower manifold of 18 bands, the Ba vibrations 
lie mostly within 100 cm~^. The in-plane Ti character 
extends throughout this lower manifold, whereas the out- 
of-plane Ti modes are mainly confined to the two narrow 
bands within 165 and 225 cm~^. Most of the Sb weight 
lies below 145 cm~^, and the out-of-plane O modes lie 
between 255 and 295 cm^^. 

The strength of the interaction between electrons and 
phonons is generally given in terms of the Eliashberg 



N{Ef) ^ ^ 



(1) 



where N(Ep) is the electronic density of states at the 
Fermi energy, eJJ is the electronic energy at wavevector 
k and band index ri, uji^q is the energy of a phonon with 
wavevector q and branch index and |5k'k+ql^ 
trix element for an electron in the state |nk) scattering to 
|mk -I- q) through a phonon cj^q. The calculated Eliash- 
berg spectral function for BaTi2Sb20 is shown in Fig. |51 
The spectral weight is spread through all the phonon fre- 
quencies, but it is specially enhanced at low frequencies 
below 90 cm~^ and between 100 and 160 cm^^. I find 
that substantial part of the spectral weight comes from 
modes that involve in-plane Ti vibrations. This is rea- 
sonable as the states near the Fermi level are dominated 
by Ti d character. Furthermore, the inter-layer Ti-Ti dis- 
tance is much larger than the intra-layer Ti-Ti distance, 
and the coupling to the out-of-plane Ti modes is much 
smaller, as one might expect. 

The q dependent total electron-phonon coupling Aq 
(= duia'^F{uj, q}/uj) is plotted for the kz ^ plane of 
the Brillouin zone in Fig. [T] There are contributions to 
the electron-phonon coupling from throughout the Bril- 
louin zone, but the magnitude of the electron-phonon 
coupling away from the zone center is much larger. In 
particular, the electron-phonon coupling is peaked near 
the wavevector (7r/a,7r/a). Again, this is expected from 
the nesting between the Fermi sheets near the zone center 
and zone corners. Furthermore, since the disconnected 
Fermi sheets have different orbital character weights, the 
electron-phonon coupling on different Fermi sheets are 
likely to be different. Hence, this material is very likely 
to host a multiband superconductivity. 

The total electron-phonon coupling constant is given 
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FIG. 7: Calculated q dependent total electron-phonon cou- 
pling Aq shown for the fcz = plane of the Brillouin zone. 
The area of the circles is proportional to the magnitude of 
Aq. The unit of the axes is 2-k /a. 



by Acp = En,^ \m = 2 /j, dua'^F{oj)/Lj. For 
BaTi2Sb20, I obtain Acp = 1.28 with the logarithmi- 
cally averaged frequency win = 65 cm^^. (I obtain these 
numbers by setting the contribution from the imaginary 
frequencies to zero.) One can estimate the Tc within 
a single-band picture by inserting these numbers in the 
simplified Allen-Dynes formula, 



1.2 



exp 



1.04(1 + A, 



■cp/ 



Acp-Af*(l + 0.62A, 



With a value for the Coulomb pseudopotcntial parameter 
/i* = 0.1, I obtain Tc = 9.0 K, which shows that the 
conventional electron-phonon picture readily explains the 
superconductivity in BaTi2Sb2 0. 

The electron-phonon calculations indicate that the su- 
perconductivity is of multiband strong-coupling nature 
with the predicted Tc of 9.0 K overestimating the experi- 
mentally obtained values of 5.5 K and 1.2 This dis- 
crepancy may be due to the pair breaking effects of spin 
fluctuations. The total coupling relevant for supercon- 
ductivity is A = Acp — Asf, where Agf is the contribution 
due to spin fluctuations. The presence of spin fluctu- 



ations, as shown to exist in this compound by Singh^ 
would therefore reduce the Tr- 



V. CONCLUSIONS 

In summary, I have presented the results of phonon 
dispersions and electron-phonon coupling calculations for 
BaTi2Sb20. The phonon dispersions show a weak lattice 
instability around Brillouin zone corners that gives rise to 
a CDW instability. This CDW phase is likely the cause 
of the experimentally observed transition at 54 K that 
is seen in the measurements of susceptibility, resistivity, 
and heat capacity measurcmcntsj^i^ The distortions cor- 
respond to elongation or compression of the Ti squares 
without an enclosed O such that the Ti squares with O 
rotate either clockwise or counterclockwise by a small 
amount. This results in a (^/2 X \/2 X 1) doubling of the 
unit cell. 

The electron-phonon calculations reveal the presence 
of strong electron-phonon couplings, especially to the 
in-plane Ti modes. There are contributions to the to- 
tal coupling from throughout the Brillouin zone, but the 
magnitude of the couplings are large away from the zone 
center. The electron-phonon coupling is peaked near the 
zone corners at wavevectors that would correspond to 
the nesting vectors between the Fermi sheets around the 
zone center and zone corners. The total electron-phonon 
coupling is Aep = 1.28, which gives an estimate of Tc = 
9.0 K when inserted in the Allen-Dynes formula. This 
overestimates the experimentally observed values of Tc- 
This discrepancy may be due to the pair breaking effects 
of spin fluctuations, which has been shown to exist in this 
compound 1^ 
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